Tuning spin-orbit coupling and superconductivity at the SrTiOs/LaAlOa interface: 

a magneto-transport study 
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The superconducting transition temperature, T c , of the SrTiOa/LaAlOs interface was varied by 
the electric field effect. The anisotropy of the upper critical field and the normal state magneto- 
transport were studied as a function of gate voltage. The spin-orbit coupling energy eso is extracted. 
This tunable energy scale is used to explain the strong gate dependence of the mobility and of the 
anomalous Hall signal observed, eso follows T c for the electric field range under study. 

PACS numbers: 75.70.Cn, 73.40.-c 



Interfaces between strongly correlated oxides exhibit a 
variety of physical phenomena and are currently at the 
focus of intensive scientific research. An electronic recon- 
struction occurring at the interfaces may be at the origin 
of these phenomena. [lj It has been demonstrated that 
the interface between SrTi0 3 (STO) and LaA10 3 (LAO) 
is highly conducting, having the properties of a two di- 
mensional electron (2DEG) gas. [2j At low temperatures 
the 2DEG has a superconducting ground state, whose 
critical temperature can be modified by an electric field 
effect. Q The origin of the charge carriers and the thick- 
ness of the conducting layers are still under debate. [H,I3- 
IH Recently, we have demonstrated that for carrier con- 
centrations at the range of 10 13 cm~ 2 , a large, hig hly 
anisotropic magnetoresistance (MR) is observed. [10j Its 
strong anisotropy suggests that it stems from a strong 
magnetic scattering confined to the interface. 

Here we show that both superconducting and spin- 
orbit (SO) interaction can be modified by applying a 
gate voltage. The upper critical field applied parallel to 
the interface, H c2 || , is approximately the weak coupling 
Clogston-Chandrasekhar paramagnetic limit: /io-ff C 2 — 
1.75fcsT c for high carrier concentrations. However, as the 
density of the charge carrier is reduced, H c2 || becomes as 
large as five times this limit, suggesting a rapidly increas- 
ing SO coupling. This SO coupling energy eso manifests 
itself in many of the transport properties studied. 

We use a sample with 15 unit cells of LaA103, de- 
posited by pulsed laser on an atomically flat SrTiO3(100) 
substrate. The oxygen pressure during the deposition 
was maintained at 10~ 4 Torr. The deposition was fol- 
lowed by a two-hour annealing stage at oxygen pressure 
of 0.2 Torr and a temperature of 400 C. The thickness 
of the LaAlOa layer was monitored by reflection high en- 
ergy electron diffraction. Growth procedure is similar to 
the one described elsewhere. A layer of gold was evap- 
orated at the bottom of the sample and used as a gate 
when biased to ±50 V relative to the 2DEG. Contacts 
were made using a wire bonder in a Van Der Pauw ge- 
ometry. We took extra care to ensure the absence of heat- 
ing by testing the resistance to be current-independent 



FIG. 1: Color on-line (a) Normalized sheet resistance versus 
temperature for the various gate voltages (from left to right) 
V 9 =50 V, as-grown state, 10 V, -10 V, -50 V color and shape 
code applies to all graphs in the paper, (b) Resistivity versus 
the perpendicular magnetic field . Inset: zoom on the low field 
region where superconductivity is suppressed. 



at the superconducting transition point, where current 
sensitivity is maximal. During measurements, the sam- 
ple was cooled down to a base temperature of 20 mK and 
a magnetic field of up to 18 T was applied. Samples ro- 
tation was done using a step motor with a resolution of 
0.015°/step. All transport properties reported here were 
measured in the as- grown state (AG), and while applying 
V s =50, 10, -10, -50 V. In order to accurately determine 
T c , the remnant magnetic field was minimized by oscil- 
lating the field down to zero. In addition, the sample was 
kept parallel to the field where superconductivity is less 
sensitive to it. 

FigQJi presents the normalized resistance R/R„ as a 
function of temperature for the various V g values. R„ 
is the resistance measured at zero field and at T=0.5 K, 
well above T c , for each V g . Tuning V g from +50 to -50 
V increases T c from 0.1 to 0.35 K. 

FigHb presents the resistance normalized with R n 
versus the perpendicular magnetic field at T=20 mK. 
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FIG. 2: Sheet resistance versus the magnetic field applied 
parallel to the interface and current for various charge densi- 
ties. 



A large, positive MR is observed for V g =50 V. Its mag- 
nitude gradually decreases as V g is reduced. The inset 
zooms on the superconducting low field region. 

Table 1 summarizes the superconducting properties for 
the various V g . We define H c2 j_ (T c ) as the field (tem- 
perature) at which the resistance is R n /2. The Ginzburg 
Landau coherence length £gl = \f < f ) o/27r/xo-ff c _L is ex- 
tracted from H c2 ^, here $o is the flux quantum. 

The charge carrier density n = I /Rue is inferred from 
Hall measurements in the high field regime (fJ-oH >14 
T). It is presented together with the calculated mobility 
/i = n/R n . A large variation of \i is observed, indicating a 
significant change in the scattering rate with gate voltage 
as previously noted. ll| This behavior will be discussed 
later. 

FigEl presents the sheet resistance versus the mag- 
netic field, H || , applied parallel to the interface and to 
the current. From the low field regime we extract the su- 
perconducting parallel critical field H c2 ||=H(R n /2). We 



note that while T c increases merely to 0.35 K, H c2 |i be- 
comes as large as 2.5 T [table 1]. The high H c2 || and low 
T c implies that the paramagnetic limit is exceeded. |12j| 
Since £gl is rather large, 50-200 nm, it is reasonable to 
use the weak coupling BCS approximation. Accordingly, 
we expect superconductivity to exist in fields lower than 
giJ-sH < 3.5fcsT c ; here g is the gyromagnetic ratio and 
fiB is the Bohr magneton. 

In FigJSJi ^sH c2 || is plotted against fcsT c for differ- 
ent V g values. The dashed straight line represents the 
paramagnetic limit using BCS weak coupling and g=2. 
As T c increases, the limit is exceeded by up to a fac- 
tor of ~ 5. This behavior suggests a strong SO coupling 
that relaxes the Clogston-Chandrasekhar limitations. 12 1 
For a superconducting layer thickness d <C £gl it is ex- 



pected that H C 2|| be much larger than H c2 j_. In this case 
the critical field is determined by the paramagnetic limit 
and the spin-orbit energy. From FigJSJi it appears that 
even for the lowest T c , superconductivity is quenched 
in this paramagnetic limit. We can therefore merely 
set an upper limit on the thickness of the conducting 
layer by analyzing the anisotropy of the critical field: 
d < v / 3 < &o/7tCglMoHc2||- 13j The thickness upper limit 



is presented for the various V g reaching a minimal value 
of 10 nm for n = 3 x 10 13 cm~ 2 [table 1]. Previous esti- 
mations obtained similar values. They, however, relate d 
to the actual thickness. 
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Let us describe the general behavior of R(Hii) [FigJ2] 



using the 



6.1 x 10 



13, 



curve (green triangle) 



as an example: Above the superconducting transition, 
the sheet resistance reaches a roughly field-independent 
regime (0.9-1.8 T) in which the resistance approximately 
equals R„. We define H* as the onset field where dR/dH 
becomes negative. For H>H* the resistance drops to 
R sat , the low resistance saturation regime (H||>10 T). 
Both H c2 || and H* strongly depend on V g ; they in- 
crease as V g changes from 50 to -50 V, and H* becomes 
unmeasurably high for V g =-50 V (black squares). 
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Figl3K suggests that SO coupling plays a major role in 
the system. We shall now describe the data presented 
so far using a single energy scale, eso- We relate H* to 
the breakdown field of the spin-orbit coupling energy: 
g/isH* = eso- The eso values are presented assum- 
ing g=2 [table 1]. Above this field SO scattering is sup- 
pressed and completely vanishes at the high field satura- 
tion regime where all spins are aligned. In this scenario 
the high field saturation value R sa t is the remnant impu- 
rity scattering. The SO scattering rate can be evaluated 
from the difference between R ra and K sa t . Therefore R„- 
H S at should be proportional to h/rso = e SO- 

FiglSb presents R„ — R sa t and /xbH* versus ksT c . 
These two quantities scale as predicted. Furthermore, 
the saturation resistance R sa t roughly scales with the 
number of carriers. We further use the calculation of 
R. A. Klemm et at which estimates the SO scattering 
time from H c2 || and TV for a two-dimensional film in the 
paramagnetic limit. [161] As seen [Fig[3]D, red circles], the 
general behavior is described by this (probably oversim- 
plified) model. We can therefore conclude that a single 
energy scale, tuned by V g , dominates the behavior of the 
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transport and affects superconductivity. For the carrier 
concentration range under study R SQ t <C R n - Hence, the 
main contribution to the zero field resistance comes from 
SO scattering. This explains the strong dependence of 
the mobility on gate voltage. 

We note that for V g = 50 V the carrier concentra- 
tion is lower compared to the as-grown state. This is in 
contrast to a simple capacitor-like behavior with nega- 
tive charge carriers. Furthermore, despite the decrease 
in carrier concentration, R„ unexpectedly decreases. We 
relate this peculiar behavior to the SO scattering pro- 
cesses dominating R„. We assume that static positive 
charges move to the interface when the highest positive 
voltage is applied and consequently n is reduced. This 
charge movement reduces the initial electric field at the 
interface and as a result decreases the SO scattering. 
Upon reducing the gate voltage from its maximal pos- 
itive value, electrons are drawn away from the interface 
in a reversible manner as expected, and an electric field 
is created. The initial electric field in the as-grown state 
results in enhanced SO scattering and higher resistance, 
despite the higher carrier concentration. 

The MR presented [FigQ}}, Fig|2] is highly anisotropic 
for high magnetic fields. For example, for n — 6.1 x 
10 13 cm -2 and H=18 T, the resistance vary from 0.12R„ 
to 3.5R„ (a factor of 30) for the in-plane and out-of- 
plane field configurations respectively. We study this 
anisotropy by rotating the sample around an in-plane 
axis which is perpendicular to the current, while keep- 
ing the total field constant, \H\ =18 T. FigfJJi presents 
the normalized sheet resistance versus the perpendicular 
field component. For the small angle range presented, 
the parallel field component is approximately constant 
~18 T. Moreover, for this parallel field component, the 
resistance is insensitive to small changes in Hm [Fig{2]. 
As shown, the sheet resistance rapidly increases when a 
small perpendicular component is introduced. It reaches 
R„ for Hx— 1-5 T (an angle of about 4°). This is sim- 
ilar to our previous observation at 2 K.[10( The curves 
n = 7.8, 6.1 x 10 13 cm -2 and AG (not shown) merge near 
1.5 T, while n = 4.4 x 10 13 cm -2 departs form the gen- 
eral behavior. This can be easily explained since for this 
gate voltage the resistance does not saturate up to 18 
T [FigE]. For n = 3.0 x lO^cm' 2 , H=18 T is well be- 
low H* and therefore the anisotropy cannot be observed 
at this too small a magnetic field. We relate the sud- 
den appearance of resistance with a small perpendicular 
field component to the orbital motion generated by this 
component. This motion is enhanced by the high mo- 
bility state (R~R sa t) induced by the presence of a large 
parallel magnetic field. 

FigHb presents the Hall resistivity p xy at 20 mK 
versus the magnetic field after subtracting out the lin- 
ear term obtained from a fit to the high field regime. A 
conspicuous deviation from this linear part in the Hall re- 
sistivity (anomalous Hall effect AHE), is observed. This 





(a) 




O 


p> u. 1U - 














O 


0.05- 






2^ B H c =1.75kJ c _ 






.-e- © — - 





2- 


(h) 


i 












a 

W o 1- 


A 


2n B H* 


o 




o 


Theory 


3 




□ 


AR R 


0- 









0.9 



0.01 0.02 0.03 

k B T c (meV) 



FIG. 3: (Color on-line) a. The parallel upper critical field H ii 
versus k_gT c while V 9 is varied. The dashed line represents 
the expected behavior for the paramagnetic limit, 
b. eso extracted from H* [see Fig[2], and calculated from 
the measured superconducting properties (H c2 |i and T c ),[? ] 
versus ksT c . The right axes presents AR, which scales with 
eso as expected. 



AHE persists up to 100 K with a peculiar, roughly linear 
temperature dependence (not shown) . As shown, the ap- 
plied V g varies the AHE saturation field. This variation 
is inconsistent with a simple magnetic impurity scenario 
and with a simple ferromagnetic behavior. Furthermore, 
no hysteresis is observed for all carrier concentrations. 
Here we note that the AHE saturation field roughly fol- 
lows the behavior of eso and may be related to it. 

FigfJJ: compares the data in Fig0t> in the case where 
n = 4.4 x 10 13 cto -2 (red circles), with data taken at 
constant magnetic field while rotating the sample (as in 
FigHk)- For the rotation (brown pluses), the slope of 
the AHE is much steeper. The two measurements dif- 
fer at low perpendicular fields where the mobility is an 
order of magnitude higher for the rotation. The qual- 
itative variation between the measurements is expected 
since the AHE is usually proportional to some power of 
the resistivity. The AHE has been related to a multiple 
band structure (light and heavy) [11] and to incipient dis- 
ordered magnetism induced by the field (l7^. The former 
is unlikely in view of the strong AHE slope in the case of 
the rotation experiment. For the case of light and heavy 
charge carriers, one would expect the former to domi- 
nate the resistivity and the latter to dominate the low 
field Hall. However, both the resistivity and the AHE 
are very susceptible to a small perpendicular field. 

Finally, we note that for the highest mobility state 
achieved (blue stars), quantum oscillations are observed 
in both MR and Hall measurements [FigdJa, Fig0J)]. 
These oscillations will be analyzed separately. [18j 

Spin-orbit interactions are expected to play a signif- 



4 




H„H ± (T) H H,(T) 

FIG. 4: (Color on-line) a. Normalized resistance as a function 
of perpendicular magnetic field. Data is taken while rotating 
the sample in a constant magnetic field of f8 T. b. Hall re- 
sistivity after subtraction of a linear fit for high fields (AHE). 
c. The AHE for n = 4.4 x fO 13 cm" 2 while H is scanned, and 
while rotating the sample at a constant field. 

icant role at the interface due to the inversion symme- 
try breaking. The SO interaction term is of the form: 
e so = p • er x VU where p is the carrier momentum, a is 
the Pauli spinor, and W is the electric field perpendic- 
ular to the interface in our system. [l9| The strong gate 
dependence of eso directly follows from this expression. 
As pointed out above, the initial local electric field is 
unknown. When a positive gate voltage was applied for 
the first time the number of charge carriers decreased, 
while the resistivity decreased. We believe that this is 
mainly due to a decrease in the total electric field near 
the interface, resulting in a decreased SO scattering and 
an enhanced mobility. This suggests that the initial field 
(AG state) is a consequence of the electronic reconstruc- 
tion. In this scenario, adding more LaA103 layers in- 
creases the as-grown field and hence the SO scattering. 
This explains the decrease in mobility while increasing 
the LaA10 3 thickness.^ 

It seems that most scattering processes in the normal- 
state involve a spin flip, which strongly impede the trans- 
port. In a simple metal such processes do not contribute 
much to the resistivity, yet in our system they become 
important due to the strong SO coupling. This is clear in 
the case of in-plane spin flip processes; since SO coupling 
affects in-plane spin orientations, a spin flip process re- 
sults in momentum reversal and consequently in a strong 
contribution to the resistivity. These spins are strongly 
coupled to the in-plane momenta; therefore it is impos- 
sible to align them with Hp unless the energy associated 
with the field fl^sHii exceeds eso- F° r H>H*, the spins 
gradually align along the field direction and so R„ is re- 
duced. When all spins are aligned, the resistance reaches 
the saturation value R sa t- For the out-of-plane field ori- 
entation, suppression of spin scattering is overwhelmed 
by the large positive MR. 



Tuning the gate voltage from 50 to -50 V increases the 
local electric field, resulting in an increase of eso an d 
H*. Unlike R n (the zero field resistance), R sat roughly 
scales with number of charge carriers deduced from the 
high field Hall measurements. This suggests that R, sa t 
is a result of standard impurity scattering. For Vg=-50 
V, eso is extremely large and an in-plane field of 18 T 
is not enough to align the spins and suppress the spin- 
scattering resistivity. Upon applying a small perpendicu- 
lar field component and due to the large scattering time, 
orbital motion is immediately turned on along with the 
SO coupling. This results in a full recovery of the spin- 
MR. 

In summary, we study the phase diagram of the 
SrTiOa/LaAlOa interface in the region where T c in- 
creases, while reducing the carrier concentration by vari- 
ation of gate voltage. We demonstrate the important 
effect of spin-orbit (SO) interaction on both supercon- 
ductivity and on normal-state transport. The SO cou- 
pling energy (eso) is evaluated using two independent 
transport properties, and is also in agreement with theo- 
retical model given the superconducting parameters: T c 
and the upper critical parallel field, eso follows T c for 
the electric field range studied. Our results suggest that 
SrTi03 /LaA103 interfaces may be useful for future oxide 
based devices controlling the orbital motion of electrons 
by acting on their spins. [2l| 
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